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Abstract 

A measurement is presented of the ratio of the inclusive 3-jet cross section to 
the inclusive 2-jet cross section as a function of the average transverse momen- 
tum, (pti,2)/ of the two leading jets in the event. The data sample was collected 
during 2011 at a proton-proton centre-of-mass energy of 7 TeV with the CMS de- 
tector at the LHC, corresponding to an integrated luminosity of 5.0 fb^ . The 
strong coupling constant at the scale of the Z boson mass is determined to be 
as(Mz) = 0.1148 ± 0.0014 (exp.) ± 0.0018 (PDF)+j^;gg^g (scale), by comparing the ratio 
in the range 0.42 < (pti,2) < 1.39 TeV to the predictions of perturbative QCD at 
next-to-leading order. This is the first determination of as(A^z) from measurements 
at momentum scales beyond 0.6 TeV. The predicted ratio depends only indirectly on 
the evolution of the parton distribution functions of the proton such that this mea- 
surement also serves as a test of the evolution of the strong coupling constant beyond 
0.42 TeV. No deviation from the expected behaviour is observed. 
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1 Introduction 

As a consequence of the non-Abelian nature of quantum chromodynamics (QCD), the renor- 
mahzation group equation (RGE) lU-El predicts that the strong force becomes weaker at short 
distances corresponding to large momentum transfers, a property of QCD referred to as asymp- 
totic freedom. The strength of the strong force, a.s{Q), at a given distance or momentum scale Q 
is not predicted and has to be extracted from experiment. Measurements at different Q can then 
be compared for consistency with QCD via the RGE, w^hich precisely describes the evolution of 
*s(fir)/ where the renormalization scale jir is identified with Q. By convention, the consistency 
is tested by evolving all values of a.s{Q) to the common scale }ir = Q = Mz, i.e. the precisely 
known mass of the Z boson. The current world average value is 0.$ (Mz) = 0.1184 ± 0.0007 [4J. 

Measurements of the running of as(Q) provide a stringent test of QCD. Previous collider ex- 
periments at LEP and HERA have established the validity of the RGE up to momentum trans- 
fers Q of 208 GeV [4J. A recent publication by the DO Collaboration extends this range up to 
400 GeV |5|. The determination of o-siQ) from jet cross sections as in [61 or [T'l depends directly 
on parton distribution functions (PDFs) that have been evolved from small to very high mo- 
mentum scales via the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations IHHTOl, 
which assume the validity of the RGE. This dependence on the evolution of the PDFs can be 
reduced by investigating cross section ratios. The ratio R32 of the inclusive 3-jet cross section 
to the inclusive 2-jet cross section is proportional to ix.s{Q) where Q is defined as the average 
transverse momentum of the tw^o jets leading in pj, 

Q = (PT1,2) = J^- (1) 

Many theoretical systematic uncertainties related to the choice of the renormalization and fac- 
torization scales, fir and ]if, or to nonperturbative effects are reduced in the cross section ratio. 
In addition, experimental uncertainties such as those due to the jet energy scale largely cancel 
in the measurement of R32. The uncertainty on the integrated luminosity measurement can- 
cels completely. The Compact Muon Solenoid (CMS) Collaboration has previously measured 
R32 fll], and the predictions of various Monte Carlo (MC) event generators w^ere found to be 
in general agreement with the measurement. 

This measurement is performed using a sample of multijet events, collected during 2011 by 
the CMS experiment at the Large Hadron Collider (LHC), corresponding to an integrated lu- 
minosity of 5.0 fb^ of pp collisions at a centre-of-mass energy ^^s = 7TeV. The transverse 
momentum pj and the rapidity 1/ of a jet with energy E and momentum p = [px, Py, Pz) ('where 
Pz is the momentum component along the direction of the anticlockwise proton beam) are de- 
fined as pj = vpl + Py and y = \ ln[(£ + pz)/{E — Pz)], respectively. Jets are reconstructed 
using the infrared- and collinear-safe anti-fcx clustering algorithm l|T2l[T3ll with a size parameter 
of 0.7. This measurement uses jets with pi > 150 GeV and \y\ < 2.5. 

The large number of multijet events collected over a wide range of (pti,2)/ 420 < (pti,2) < 
1390 GeV, allows ois{Q) to be determined with only a small dependence on the evolution of the 
PDFs, thus testing the validity of the RGE in an extended range of transverse momenta up to 
ITeV. 



2 3 Event selection and reconstruction 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, 13 m in length and 
6 m in diameter, providing an axial magnetic field of 3.8 T. The field volume of the solenoid 
is instrumented with various layers of particle detection systems. Charged particle trajecto- 
ries are measured by the silicon pixel and strip tracker, with full azimuthal coverage within 
\ri\ < 2.5, where the pseudorapidity t] is defined ast] = — ln[tan(0/2)], and 6 is the polar angle 
with respect to the z axis. Surrounding the trackers are a lead tungstate crystal electromag- 
netic calorimeter (ECAL) with a preshow^er detector in the endcaps, and a brass and scintillator 
hadron calorimeter (HCAL), covering the region |/y| < 3. In addition to the barrel and endcap 
detectors, CMS has extensive forward calorimetry which extends the coverage to | // 1 = 5. The 
steel flux return yoke outside the solenoid is instrumented with gas-ionization detectors used 
to identify and reconstruct muons. A more detailed description of the CMS detector can be 
found in [14 J. 

3 Event selection and reconstruction 

The CMS detector records events using a two-level trigger system consisting of a hardw^are- 
based level-1 (LI) trigger and a software-based high level trigger (HLT). In this study, single- 
jet triggers that reconstruct jets from calorimeter energy deposits at LI and HLT are used to 
select events based on three HLT -pj thresholds, 190, 240, and 370 GeV. All except the highest- 
threshold trigger were prescaled during the 2011 run. The corresponding integrated luminosity 
C for each of the three samples is shown in Table IT] The efficiency of each of the triggers is 
estimated using lower-px"threshold triggers. These three jet trigger thresholds ensure 100% 
trigger efficiency in the three jet samples for (pti,2) > 215, 269, and 409 GeV. 

Table 1: The integrated luminosity for each trigger sample. 
HLT pt threshold (GeV) 190 240 370 



C{{b'') 0.15 0.51 5.0 



Each event is required to have at least one offline-reconstructed vertex llTSi along the beam 
line that is within 24 cm of the nominal interaction point. The four-vectors of particle candi- 
dates reconstructed by the CMS global event reconstruction algorithm (also called particle-flow 
event reconstruction [16J) are used as input to the jet-clustering algorithm. The clustering is 
performed by the FastJet package [13] using four-momentum summation. The global event 
reconstruction algorithm reconstructs and identifies each particle with an optimized combina- 
tion of subdetector information. The energy of photons is obtained directly from the ECAL 
measurements after being corrected for zero-suppression effects. The energy of electrons is 
determined from a combination of the track momentum at the main interaction vertex, the cor- 
responding ECAL cluster energy, and the energy sum of all bremsstrahlung photons attached 
to the track. The energy of muons is derived from the corresponding track momentum. The 
energy of charged hadrons is determined from a combination of the track momentum and 
the corresponding ECAL and HCAL energies, corrected for zero-suppression effects, and cali- 
brated for the nonlinear response of the calorimeters. Finally, the energy of neutral hadrons is 
obtained from the corresponding calibrated ECAL and HCAL energies. 

Jet energy corrections [17J are derived using simulated events, generated by pythia 6.4.22 flSl 
and processed through the CMS detector simulation based on Geant4 [19J, and in situ mea- 



surements with dijet, photon+jet, and Z+jet events. An offset correction is applied to take into 
account the extra energy clustered into jets from additional proton-proton interactions within 
the same or neighbouring bunch crossings (in-time and out-of-time pileup) [17J. Pileup effects 
are important only for low-px jets and become negligible for jets with pj > 200 GeV. The cur- 
rent measurement is therefore largely insensitive to pileup effects. The jet energy corrections, 
which depend on the // and pj of the jet, are applied to the jet four-momentum vector as a 
multiplicative factor. The multiplicative factor is in general smaller than 1.2, approximately 
uniform in rj, with typical values of 1.1 for jets having pj = 100 GeV and decreasing to 1.0 for 
higher values of pj. 

To suppress nonphysical jets, i.e. jets resulting from noise in the ECAL and/or HCAL calorime- 
ters, tight identification criteria are applied: each jet should contain at least two particles, one 
of which is a charged hadron, and the jet energy fraction carried by neutral hadrons, photons, 
muons, and electrons should be less than 90%. These criteria have an efficiency greater than 
99% for genuine jets. 

The selection of multijet events requires two or more jets w^ith transverse momentum greater 
than 150 GeV and \y\ < 2.5. The final sample is extracted by rejecting events if either or both of 
the leading jets in pj have \y\ > 2.5. 

4 Measurement of 1^32 and comparison with theoretical predic- 
tions 

The measured ratio R32 as a function of (pti,2) is the ratio of the number of selected inclusive 
3-jet events to the number of selected inclusive 2-jet events in each (^11,2) bin. The ratio R32 
is corrected for detector smearing effects and unfolded to stable-particle level. The unfolding 
method is the iterative Bayesian method [20 1, as implemented in the RooUnfold software 
package 121]. Unfolding uses a response matrix that maps the true distribution onto the mea- 
sured one. The response matrix is derived from a simulation, which uses as input the true 
R32 distribution from PYTHIA6 tune Z2 and introduces the smearing effects by taking into ac- 
count the (pti,2) resolution IllTl . After unfolding R32 to stable-particle level, the final statistical 
uncertainties include the correlation among the various (^11,2) bins. 

Two main sources of systematic uncertainties on R32 are considered. The first is due to the jet 
energy scale (JES) and the second due to the unfolding. 

The JES uncertainty has been estimated to be 2.0-2.5% for particle-flow jets p2| . depending on 
the jet Pj and rj. All mutually uncorrelated JES uncertainty sources are considered following 
the procedure described in Ref. |23|. The total systematic uncertainty on R32 due to the JES 
uncertainty is 1.2%. 

The unfolding method takes into account three different mutually uncorrelated uncertainty 
sources. The first arises from insufficient knowledge of the simulated inclusive 3-jet and 2-jet 
(^11,2) spectra, which are used to construct the simulated ratio R32 used in the unfolding. The 
uncertainty is estimated by varying the 3-jet and 2-jet spectra slopes by ±10%. This is a conser- 
vative estimate and is motivated by the observed difference in the 3-jet and 2-jet spectra slopes 
between simulations using the event generators PYTHIA6 tune Z2 and HERWIG++ ||24| with its 
default tune in version 2.3. Simulations of 3-jet and 2-jet spectra using the MadGraph I^25ll26| 
event generator are in agreement w^ith those of PYTHIA6 . The second uncertainty arises from 
the insufficient knowledge of the (pti,2) resolution and is estimated by varying it by ±10%. 
This variation is motivated by the observed difference between data and simulation in the jet 



4 4 Measurement of R32 and comparison with theoretical predictions 

energy resolution IllTl . Finally, the third uncertainty arises from non-Gaussian components 
in the (^11,2) resolution and is estimated by adding non-Gaussian tails to the simulation. The 
overall systematic uncertainty on R32 due to unfolding is less than 1%. A potential bias originat- 
ing from the unfolding technique is studied by comparing the unfolding result of the Bayesian 
method with that of the singular- value decomposition (SVD) method l|27| . The bias is found to 
be negligible. 

The theoretical calculation of the ratio R32 is based on next-to-leading-order (NLO) perturba- 
tive QCD (pQCD) calculations multiplied by a nonperturbative factor, which corrects for mul- 
tiparton interactions (MPI) and hadronization effects. The NLO calculations assume Nf = 5 
massless quark flavours and are based on the parton-level generator NLOJET++ Il28ll29ll . A 
test exaggerating the impact of the top quark by imposing Nf = 6 massless flavours revealed 
only negligible changes at the percent level in the ratio. The computations with NLOJET++ 
are performed w^ithin the FASTNLO framew^ork |30l |3TJ using the following four PDF sets: 
NNPDF2.1 13211331, ABMll iH, MSTW2008 ESIIMI, and CTIO ISZIIMI- In each case the PDF 
version employing next-to-next-to-leading order (NNLO) evolution code is chosen, and for 
comparisons the respective default values of 0.$ (Mz), which are 0.119, 0.1134, 0.1171, and 0.118, 
are used. NNPDF2.1, MSTW2008, and CTIO utilize a variable-flavour-number scheme with the 
maximal number of flavours Nf max equal to 6, 5, and 5 respectively, while the ABMll PDF set 
was developed in a fixed-flavour-number scheme with Nf = 5. The renormalization and fac- 
torization scales are set to the average transverse momentum (^11,2)- 

The nonperturbative effects are estimated using the PYTHIA6 tune Z2 and HERWIG++ tune 2.3 
event generators. The chosen MC models feature different descriptions of the phenomena and 
are representative of the possible values of the nonperturbative corrections. The nonperturba- 
tive correction (NPC) factor is defined as the ratio of R32 predicted with the nominal generator 
settings to that obtained w^ith the MPI and hadronization switched off. This factor is calcu- 
lated considering the average of the two MC generators and has typical values of ?»1.02 for 
(Pti,!) = 250 GeV, decreasing to 1.0 for higher (pti,2)- The uncertainty is considered to be half 
of the difference between the NPC values obtained using the tw^o MC generators and amounts 
to ?»0.1%, leading to a negligible influence on the final result. 

Finally, uncertainties due to the renormalization and factorization scales are evaluated by vary- 
ing from the default choice of jiy = }if = (^11,2) between (pxi,2)/2 and 2(pxi,2)/ simultane- 
ously in the numerator and denominator of the ratio R32, in the following six combinations: 

{^r/{p^l,2),^f/{p^l,2)) = (1/2,1/2), (1/2,1), (1,1/2), (1,2), (2,1), and (2,2). 

Figure [1] presents the measured ratio R32 together with NLO predictions using the NNPDF2.1 
(top left), the ABMll (top right), the MSTW2008 (bottom left), and the CTIO (bottom right) 
NNLO PDF sets. The upper panel of each plot shows the ratio R32 (solid circles) together with 
the NLO prediction (solid line) corrected for nonperturbative effects, the scale uncertainty, and 
the PDF uncertainty. At the bottom of each plot, the ratio of data over theory is shown together 
with bands representing the scale (dotted lines) and PDF uncertainties (solid lines). The error 
bars in the figure correspond to the total uncertainty, for which the statistical and systematic 
uncertainties are added in quadrature. For each PDF set the respective default value of as(Mz) 
is used in this comparison as indicated. 

The measured ratio rises with increasing (^11,2) as the phase space opens up for the production 
of a third jet, reaching a plateau value for 600 < (pti,2) < 1000 GeV. At higher (^11,2) R32 
decreases again because of the running of a.s, smaller parton luminosities, and because 3-jet 
configurations reach kinematic limits earlier than dijet events. 
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Figure 1: Measurement of R32 and NLO predictions using the NNPDF2.1 (top left), the ABMll 
(top right), the MSTW2008 (bottom left), and the CTIO (bottom right) NNLO PDF sets. In the 
upper panel of each plot, the ratio R32 (solid circles) together with the NLO prediction (solid 
line) corrected for nonperturbative effects (NPC), the scale uncertainty, and the PDF uncer- 
tainty are shown. The bottom panels show the ratio of data to the theoretical predictions, to- 
gether with bands representing the scale (dotted lines) and PDF (solid lines) uncertainties. The 
error bars correspond to the total uncertainty. For each PDF set the respective default value of 
(x.s{Mz) is used as indicated. 
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Scale uncertainties have a very similar behaviour for all PDF sets and dominate the region up 
to (pti,2) ~ 400 GeV. A comparison to jets with a size parameter oi R = 0.5 reveals consistent 
results but with larger scale uncertainties. 

The PDF uncertainties are different for each individual PDF set. The CTIO set exhibits the 
largest PDF uncertainties, which are of the order of 2% at (^11,2) = 400 GeV, increasing to 2.5% 
in the 1 TeV region. For the NNPDF2.1 set PDF uncertainties are of the order of 1.5% at 400 GeV, 
increasing to 2.3% at 1 TeV. Finally, for the MSTW2008 and ABMll sets PDF uncertainties are 
of the order of 1% throughout the range of this measurement. 

The comparison of data with the predictions of pQCD in Fig. IT] demonstrates that the NLO 
calculations using the NNPDF2.1, MSTW2008, and CTIO PDF sets are in agreement with the 
measured ratio R32 throughout the range of this measurement. The NLO result employing the 
ABMll PDF set underestimates R32, especially for (pti,2) < 600 GeV. 

5 Determination of ocsiMz) 

The measurement of the ratio R32 is used for the determination of the strong coupling constant 
as(Mz). Figure |2] shows the predictions using the NNPDF2.1 (top left), the ABMll (top right), 
the MSTW2008 (bottom left), and the CTIO (bottom right) NNLO PDF sets for a series of values 
of ocs{Mz), together with the measured R32. The as(Mz) value is varied in the range 0.106- 
0.124, 0.104-0.120, 0.107-0.127, and 0.110-0.130 in steps of 0.001 for the NNPDF2.1, ABMll, 
MSTW2008, and CTIO PDF sets, respectively. 

From Fig. |2]one observes that the sensitivity of the ratio R32 to variations of the strong coupling 
by Aocs (Mz) = ±0.001 is different for each of the four PDF sets. This translates into differences 
in the experimental uncertainty in the value of as(Mz) obtained for each PDF set. 

The value of as(Mz) is determined by minimizing the x^ between the experimental measure- 
ment and the theoretical predictions. The x^ is defined as 

X^ = M^C-^M, (2) 

where M is the vector of the differences between the data (R32) and the theoretical values (T32) 
in each bin i, 

M' = 42 - Tk, (3) 

and C is the covariance matrix including all experimental (statistical, JES, and unfolding) un- 
certainties. C is defined as 

C = COV^'^* + Y, COV^^S Sources ^ ^ (^Q^Unfolding Sources^ (4) 

where Cov is the statistical covariance matrix that accounts for the correlations due to un- 
folding, and CovJ^s Sources^ (^Q^unfoiding Sources ^^^ ^^^ covariance matrices that account for the 

JES and unfolding systematic uncertainty sources, respectively. Each systematic uncertainty 
source for the JES and unfolding is treated as 100% correlated across the (^11,2) bins. 

To avoid the region w^ith large scale uncertainties close to the minimal jet px requirements, 
visible in Fig.[ll as(Mz) is extracted only for (^11,2) > 420 GeV. The central result is obtained 
by minimizing the x^ (Eq. 01) with respect to dcs{Mz) for the NNPDF2.1 PDF set, which is the 
only one that permits the propagation of PDF uncertainties to the fits for each value of ag (Mz). 
The experimental uncertainties are obtained from the as(Mz) values for which x^ is increased 
by one with respect to the minimum value. The result of a fit to the region of 420-1390 GeV is 

as (Mz) = 0.1148 ± 0.0014 (exp.), (5) 
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Figure 2: The NLO predictions using the NNPDF2.1 (top left), the ABMll (top right), the 
MSTW2008 (bottom left), and the CTIO (bottom right) NNLO PDF sets for a series of values 
of as(Mz), together with the measured R32. The as(Mz) value is varied in the range 0.106- 
0.124, 0.104-120, 0.107-0.127, and 0.110-0.130 in steps of 0.001 for the NNPDF2.1, ABMll, 
MSTW2008, and CTIO PDF sets, respectively. 
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with X.'^/^doi = 22.0/20 at minimum. The experimental uncertainty contains the statistical, 
JES, and unfolding sources (Eq. Q), with the JES uncertainty being the dominant one. 



The contribution of PDFs to the uncertainty of the measurement is evaluated by repeating the 
fit for each of the 100 PDF replicas of the NNPDF2.1 set at the relevant value for as(Mz). In 
this way 100 determinations of as(Mz) are obtained, whose distribution corresponds to the 
propagation of the underlying probability density from the PDFs to the fitted strong coupling. 
The PDF uncertainty of the measurement is then computed as the standard deviation of this 
distribution. A more detailed description of the method can be found in Ref. \^]. 

The uncertainties due to the renormalization and factorization scales are treated separately by 
varying the default choice of ]ir = Hf = (Pti,2) between (pxi,2)/2 and 2(pxi,2) in six combi- 
nations as explained in Section H The x^ minimization with respect to as(Mz) is repeated for 
these six combinations. The contribution from the }ir, Hf scale variations to the uncertainty in 
the measurement is evaluated by considering the differences between the NNPDF2.1 ocg^Mz) 
central value and the highest and lowest values found in these six scale combinations. Out 
of all scale combinations the lowest a.s{Mz) value corresponds to the default scale choice of 
}ir = Hf = {pii,i) and the highest to the scale choice of f/,- = jif = (pTi,2)/2. The frequent 
observation of asymmetric scale uncertainties with larger downward uncertainties in the case 
of NLO cross sections is transformed into a purely upward uncertainty for the ratio, as can be 
seen in Tabled 

Table 2: The values of as(Mz) at the central scale and for the six scale factor combinations. 



¥r/{pjl,l) 


m/{pii,2) 


as(Mz)± (exp.) 


;tVNdof 


1 


1 


0.1148 ±0.0014 


22.0/20 


1/2 


1/2 


0.1198 ±0.0021 


30.6/20 


1/2 


1 


0.1149 ±0.0014 


22.2/20 


1 


1/2 


0.1149 ±0.0014 


22.2/20 


1 


2 


0.1150 ±0.0015 


21.9/20 


2 


1 


0.1159 ±0.0014 


20.7/20 


2 


2 


0.1172 ±0.0018 


21.3/20 



The final result is 

as (Mz) = 0.1148 ± 0.0014 (exp.) ± 0.0018 (PDF)+g;gg^{^ (scale), (6) 

in agreement with the world average value of dg (Mz) = 0.1184 ± 0.0007 [4], with the Tevatron 
results [5,^.40], and a recent result obtained with LHC data [7]. 

The determination of as(Mz), which is based on the NNPDF2.1 PDF set, is also in agreement 
with the results obtained using the MSTW2008 or CTIO PDF sets 

MSTW2008: as(Mz) =0.1141 ±0.0022 (exp.), 
CTIO: ccs (Mz) = 0.1135 ± 0.0019 (exp.), 

with x^/Ndof = 20.6/20 and 21.1/20, respectively If PDF sets with NLO evolution are used 
instead the impact on the results of the fits to the ratio observable R32 is negligible. This is 
in contrast to fits to cross sections, where NNLO PDF sets usually lead to smaller values of 



ocs{Mz) than NLO ones, and confirms the reduced dependence of R32 on details of the PDF 
evolution. 

In the case of the ABMll PDFs the series in values of 0.$ (Mz) ends at 0.120, which is insufficient 
for a derivation of the complete shape of the x^ curve at minimum such that a fit value for 
a.s{Mz) including uncertainties can only be extrapolated to give around as(Mz) = 0.1214 ± 
0.0020 (exp.) with x^/^doi = 20.6/20. For the ABMll PDF set at NLO with a default value 
of as(Mz) = 0.118, the series in iXs{Mz) values ends at 0.130 such that a fit can be performed 
which yields as(Mz) = 0.1214 ± 0.0018 (exp.), consistent with the extrapolation above. The 
fit exhibits, however, a somewhat larger value of X^/^doi = 28.5/20 compared to the other 
results. 

It is observed that with ABMll PDFs a higher value of ocs{Mz) is preferred. This is in accord 
w^ith the fact that the ABMll gluon density in the phase space relevant for this analysis is 
significantly smaller than that of all other PDF sets. Thus, the fit favors a larger as(Mz) value 
to compensate for this effect. In summary, the ABMll PDF set does not describe the data as 
well as the alternative PDF sets, as shown in Fig. Il| which leads to an inferior fit quality and a 
less consistent result for the strong coupling. 

To investigate the running of the strong coupling constant in more detail, the fitted region of 
420-1390 GeV is split into three bins of (^11,2) and the fitting procedure is repeated in each of 
these bins. The three separate extractions of as(Mz) are reported in Table Is] The experimental 
uncertainties in the three obtained values are correlated. These as {Mz) determinations are then 
evolved back to the corresponding values (x.s{Q) using the 3-loop solution to the RGE from the 
NNPDF2.1 set. For each fit region the momentum scale Q, defined to be the cross-section- 
weighted average of (^11,2) from the inclusive dijet calculation at NLO with NLOJET++, is 
computed to be Q = 474, 664, and 896 GeV, respectively. These values, derived again with the 
fastNLO framew^ork, are identical within about 1 GeV for different PDFs and vary at most by 
a few GeV when using inclusive 3-jet events. 

To emphasize that theoretical uncertainties limit the achievable precision. Table [4] presents the 
decomposition of the total uncertainty for the three bins in (pti,2) into the experimental, PDF, 
and scale components. 

Figure [3] presents the strong coupling as{Q) (solid line) and its total uncertainty (band) as 
evolved from the CMS determination, as{Mz) = 0.1148+o;oo23/ using the 3-loop solution to 
the RGE from NNPDF2.1, as before. The extractions of ci.s{Q) in three separate ranges of Q as 
presented in Table |3] are also shown. In the same figure the values of as at lower scales deter- 
mined by the HI 1,41. ,42 J , ZEUS [43 1, and DO [5, 40 1 Collaborations are shown for comparison. 
The results on «§ reported here are consistent with the energy dependence predicted by the 
RGE and extend the HI, ZEUS, and DO results to the region 420 < (pTi,2) < 1390 GeV, in 
which the RGE is tested for the first time. 

Table 3: The separate determinations of ag in bins of (pti,2)- 



(Pti,2) range 
(GeV) 


Q 
(GeV) 


as(Mz) 


'^siQ) 


No. of data 
points 


^VNdof 


420-600 
600-800 
800-1390 


474 
664 
896 


Q-l -147 +0.0061 
U.ii^/ _o.0021 

n 1 1 a9 +0.0050 

U.110^_QQg32 

1 1 70 +00058 


0936 +°-°°'*° 
0894 +°-°°3i 

U.UOV^ -0.0019 

0889 +0-0033 
U.UOO^^ -0.0018 


6 

5 
10 


4.4/5 
5.9/4 
5.7/9 
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Figure 3: The strong coupling o-siQ) (solid line) and its total uncertainty (band) evolved from 
the CMS determination ag (Mz) = 0.1148tooo23 using a 3-loop solution to the RGE as a function 
of the momentum transfer Q = (pti,2)- The extractions of ocg^Q) in three separate ranges of Q 
as presented in Table |3] are shown together with results from the HI [|4Tl 142, 1 , ZEUS [|43,| , and 
DO IHUlOl experiments at the HERA and Tevatron colliders. 
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Table 4: Uncertainty composition for a.s{Mz) from the determination of as in bins of (^11,2)- 



(Pti,2) range 
(GeV) 


Q 
(GeV) 


as(Mz) 


exp. 


PDF 


scale 


420-600 
600-800 
800-1390 


474 
664 
896 


0.1147 
0.1132 
0.1170 


±0.0015 
±0.0018 
±0.0024 


±0.0015 
±0.0025 
±0.0021 


+0.0057 
-0.0000 
+0.0039 
-0.0000 
+0.0048 
-0.0003 



6 Summary 

The ratio R32 of the inclusive 3-jet cross section to the inclusive 2-jet cross section, for jets with 
Pt > 150 GeV and \y\ < 2.5, has been measured in the range 250 < (pti,2) < 1390 GeV for 
proton-proton collisions at a centre-of-mass energy of 7 TeV. The results have been compared 
with predictions of QCD at NLO obtained with various PDF sets. The NLO calculations using 
the NNPDF2.1, MSTW2008, and CTIO NNLO PDF sets are in agreement with the measured 
ratio R32 throughout the range of {pji,!} studied. However, calculations using the ABMll PDF 
sets underestimate R32 for (^11,2) < 600 GeV. 

Measurements of R32 over the range 420 < (pti,2) < 1390 GeV have been used to determine 
the strong coupling constant 0.$ at the scale of the Z boson mass. The final result is 

as(Mz) = 0.1148 ± 0.0014 (exp.) ± 0.0018 (PDF)+g:gg^g (scale) = 0.1148 taool. 

where experimental, PDF, and scale uncertainties have been added quadratically to give the to- 
tal uncertainty. The result is in agreement with the world average value of ag (Mz) = 0.1184 ± 
0.0007 IJH and represents the first determination of the strong coupling constant from measure- 
ments at scales Q of the order of 1 TeV. The dominant uncertainties are of theoretical origin 
and limit the currently achievable precision. The predicted ratio depends only indirectly on 
the evolution of the parton distribution functions of the proton and consequently this measure- 
ment also serves as a test of the evolution of (x.s{Q) beyond 0.42 TeV. No deviation from the 
expected behaviour is observed. 
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